suppressPackageStartupMessages(library(tidyverse))
devtools::load_all('~/Google Drive/My Drive/Scripts/R_packages/myUtilities/')
## ℹ Loading myUtilities

Settings

data_dir <- '/Volumes/Mitsu_NGS_3/METTL2A/'

wd <- "~/Google Drive/My Drive/Analysis/METTL2A/"
setwd(wd)

figdir <- paste0(wd, 'Figures/DRS_m3C_RNAs/Isoforms/Mapping_m3Csites_isoforms_ecdf_compare/')
tabledir <- paste0(wd, 'Tables/DRS_m3C_RNAs/Isoforms/')

theme_set(
  theme_classic(base_size = 7) +
    theme(legend.position = 'bottom')
)

Functions

read_signal_parameter_files <- function(
    .glob, 
    .dir = '/Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/', 
    ...) {
  
  fs::dir_ls(.dir, glob = .glob) |> 
    read_tsv(
      col_names = c('ref_pos', 'ref_kmer', 'dwell_time', 'intensity'), 
      col_types = 'dcdd', 
      id = 'filename', ...
    ) |> 
    separate(filename, sep = '_',
             into = c(NA, NA, NA, NA, NA, 'sample_num', 'type', 'si', NA, 'transcript_id')) |> 
    mutate(transcript_id = str_remove(transcript_id, '.tsv'))
  
}

remove_outlier_from_eachsample <- function(df, outlier_thresh = .05) {
  df |>
    group_by(transcript_id, ref_pos, sample_num) |> 
    filter(intensity > quantile(intensity, probs = outlier_thresh / 2), na.rm = TRUE)  |> 
    filter(intensity < quantile(intensity, probs = 1 - outlier_thresh / 2), na.rm = TRUE) |> 
    ungroup()
}

Read data

tsv_files <- 
  fs::dir_ls(
    '/Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/', 
    glob = '*.tsv') |> 
  fs::path_filter(
    regexp = 'ENST00000368716.9|ENST00000354332.8|ESPRESSO:chr1:1668:8'
  )
tsv_files
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_7_siMETTL2A_I_N1_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_7_siMETTL2A_I_N1_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_8_siMETTL2A_I_N2_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_8_siMETTL2A_I_N2_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_9_siMETTL2A_I_N3_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_9_siMETTL2A_I_N3_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_4_siMETTL2A_G_N1_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_4_siMETTL2A_G_N1_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_5_siMETTL2A_G_N2_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_5_siMETTL2A_G_N2_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_6_siMETTL2A_G_N3_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_6_siMETTL2A_G_N3_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_1_Cont_D_N1_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_1_Cont_D_N1_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_2_Cont_D_N2_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_2_Cont_D_N2_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_3_Cont_D_N3_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_3_Cont_D_N3_ENST00000368716.9.tsv
parameters <- 
  tsv_files |> 
  map(read_signal_parameter_files) |> 
  reduce(bind_rows)
parameters
## # A tibble: 20,090,588 × 8
##    sample_num type     si    transcript_id ref_pos ref_kmer dwell_time intensity
##    <chr>      <chr>    <chr> <chr>           <dbl> <chr>         <dbl>     <dbl>
##  1 7          siMETTL… I     ENST00000354…       7 CTCTA       0.149        77.9
##  2 7          siMETTL… I     ENST00000354…       8 TCTAC       0.0113       84.7
##  3 7          siMETTL… I     ENST00000354…       9 CTACA       0.0150       83.5
##  4 7          siMETTL… I     ENST00000354…      10 TACAA       0.00633      81.8
##  5 7          siMETTL… I     ENST00000354…      11 ACAAC       0.00267      83.5
##  6 7          siMETTL… I     ENST00000354…      12 CAACC       0.0120       89.2
##  7 7          siMETTL… I     ENST00000354…      13 AACCC       0.0283       83.8
##  8 7          siMETTL… I     ENST00000354…      14 ACCCT       0.00733      67.9
##  9 7          siMETTL… I     ENST00000354…      15 CCCTC       0.00500      73.6
## 10 7          siMETTL… I     ENST00000354…      16 CCTCT       0.00233      78.7
## # ℹ 20,090,578 more rows

ECDF

S100A4_TCCTG <- 
  parameters |> 
  filter(ref_kmer == 'TCCTG') |> 
  filter(
    (transcript_id == 'ENST00000368716.9' & ref_pos == 55) |
      (transcript_id == 'ENST00000354332.8' & ref_pos == 100)
  )
S100A4_TCCTG |> 
  export_tsv(outdir = tabledir, compression = 'gz')
## 
## Exported to: ~/Google Drive/My Drive/Analysis/METTL2A/Tables/DRS_m3C_RNAs/Isoforms/S100A4_TCCTG_2024-07-29.tsv.gz
## # A tibble: 39,628 × 8
##    sample_num type     si    transcript_id ref_pos ref_kmer dwell_time intensity
##    <chr>      <chr>    <chr> <chr>           <dbl> <chr>         <dbl>     <dbl>
##  1 7          siMETTL… I     ENST00000354…     100 TCCTG       0.00467      78.7
##  2 7          siMETTL… I     ENST00000354…     100 TCCTG       0.00267      73.7
##  3 7          siMETTL… I     ENST00000354…     100 TCCTG       0.0200       72.8
##  4 7          siMETTL… I     ENST00000354…     100 TCCTG       0.00300      73.1
##  5 7          siMETTL… I     ENST00000354…     100 TCCTG       0.00967      74.8
##  6 7          siMETTL… I     ENST00000354…     100 TCCTG       0.0110       74.1
##  7 7          siMETTL… I     ENST00000354…     100 TCCTG       0.00633      76.3
##  8 7          siMETTL… I     ENST00000354…     100 TCCTG       0.0110       77.1
##  9 7          siMETTL… I     ENST00000354…     100 TCCTG       0.0163       71.9
## 10 7          siMETTL… I     ENST00000354…     100 TCCTG       0.0100       74.8
## # ℹ 39,618 more rows
S100A4_TCCTG_intensity_ecdf <- 
  parameters |> 
  filter(ref_kmer == 'TCCTG') |> 
  filter(
    (transcript_id == 'ENST00000368716.9' & ref_pos == 55) |
      (transcript_id == 'ENST00000354332.8' & ref_pos == 100)
  ) |>
  remove_outlier_from_eachsample() |> 
  ggplot(aes(x = intensity, colour = si, linetype = transcript_id, fill = ref_kmer)) +
  scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
  stat_ecdf(lwd = 1.1) 
S100A4_TCCTG_intensity_ecdf |> 
  ggsave_multiple_formats(
    outdir = figdir, width = 6, height = 4, fontsize = 7
  )

S100A4_ATCCT <- 
  parameters |> 
  filter(ref_kmer == 'ATCCT') |> 
  filter(
    (transcript_id == 'ENST00000368716.9' & ref_pos == 54) |
      (transcript_id == 'ENST00000354332.8' & ref_pos == 99)
  )
S100A4_ATCCT |> 
  export_tsv(outdir = tabledir, compression = 'gz')
## 
## Exported to: ~/Google Drive/My Drive/Analysis/METTL2A/Tables/DRS_m3C_RNAs/Isoforms/S100A4_ATCCT_2024-07-29.tsv.gz
## # A tibble: 39,634 × 8
##    sample_num type     si    transcript_id ref_pos ref_kmer dwell_time intensity
##    <chr>      <chr>    <chr> <chr>           <dbl> <chr>         <dbl>     <dbl>
##  1 7          siMETTL… I     ENST00000354…      99 ATCCT       0.0307       74.6
##  2 7          siMETTL… I     ENST00000354…      99 ATCCT       0.00767      71.1
##  3 7          siMETTL… I     ENST00000354…      99 ATCCT       0.00733      72.3
##  4 7          siMETTL… I     ENST00000354…      99 ATCCT       0.0123       72.5
##  5 7          siMETTL… I     ENST00000354…      99 ATCCT       0.00700      70.9
##  6 7          siMETTL… I     ENST00000354…      99 ATCCT       0.0107       73.9
##  7 7          siMETTL… I     ENST00000354…      99 ATCCT       0.00300      73.3
##  8 7          siMETTL… I     ENST00000354…      99 ATCCT       0.00433      76.5
##  9 7          siMETTL… I     ENST00000354…      99 ATCCT       0.00267      71.9
## 10 7          siMETTL… I     ENST00000354…      99 ATCCT       0.00267      72.0
## # ℹ 39,624 more rows
S100A4_ATCCT_intensity_ecdf <- 
  parameters |> 
  filter(ref_kmer == 'ATCCT') |> 
  filter(
    (transcript_id == 'ENST00000368716.9' & ref_pos == 54) |
      (transcript_id == 'ENST00000354332.8' & ref_pos == 99)
  ) |>
  remove_outlier_from_eachsample() |> 
  ggplot(aes(x = intensity, colour = si, linetype = transcript_id, fill = ref_kmer)) +
  scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
  stat_ecdf(lwd = 1.1) 
S100A4_ATCCT_intensity_ecdf |> 
  ggsave_multiple_formats(
    outdir = figdir, width = 6, height = 4, fontsize = 7
  )

Density

S100A4_TCCTG_intensity_density <- 
  parameters |> 
  filter(ref_kmer == 'TCCTG') |> 
  filter(
    (transcript_id == 'ENST00000368716.9' & ref_pos == 55) |
      (transcript_id == 'ENST00000354332.8' & ref_pos == 100)
  ) |>
  remove_outlier_from_eachsample(outlier_thresh = .1) |> 
  ggplot(aes(x = intensity, fill = si, linetype = transcript_id, colour = si)) +
  scale_fill_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
  scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
  geom_density(alpha = 1/2) +
  facet_wrap( ~ transcript_id, ncol = 1)
S100A4_TCCTG_intensity_density |> 
  ggsave_multiple_formats(
    outdir = figdir, width = 6, height = 4, fontsize = 7
  )
## Warning in grid.Call.graphics(C_polygon, x$x, x$y, index): semi-transparency is
## not supported on this device: reported only once per page

parameters |> 
  filter(ref_kmer == 'TCCTG') |> 
  filter(
    #(transcript_id == 'ENST00000368716.9' & ref_pos == 55) 
      (transcript_id == 'ENST00000354332.8' & ref_pos == 100)
  ) |>
  remove_outlier_from_eachsample() |> 
  ggplot(aes(
    x = intensity, colour = si, linetype = sample_num, fill = ref_kmer)) +
  scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
  stat_ecdf(lwd = 1.1) 

Sesison info

sessioninfo::session_info()
## ─ Session info ───────────────────────────────────────────────────────────────
##  setting  value
##  version  R version 4.3.1 (2023-06-16)
##  os       macOS Sonoma 14.5
##  system   aarch64, darwin20
##  ui       X11
##  language (EN)
##  collate  en_US.UTF-8
##  ctype    en_US.UTF-8
##  tz       Asia/Tokyo
##  date     2024-07-29
##  pandoc   3.1.1 @ /Applications/RStudio.app/Contents/Resources/app/quarto/bin/tools/ (via rmarkdown)
## 
## ─ Packages ───────────────────────────────────────────────────────────────────
##  ! package      * version    date (UTC) lib source
##    abind          1.4-5      2016-07-21 [1] CRAN (R 4.3.0)
##    backports      1.4.1      2021-12-13 [1] CRAN (R 4.3.0)
##    bit            4.0.5      2022-11-15 [1] CRAN (R 4.3.0)
##    bit64          4.0.5      2020-08-30 [1] CRAN (R 4.3.0)
##    broom          1.0.5      2023-06-09 [1] CRAN (R 4.3.0)
##    bslib          0.5.1      2023-08-11 [1] CRAN (R 4.3.0)
##    cachem         1.0.8      2023-05-01 [1] CRAN (R 4.3.0)
##    callr          3.7.3      2022-11-02 [1] CRAN (R 4.3.0)
##    car            3.1-2      2023-03-30 [1] CRAN (R 4.3.0)
##    carData        3.0-5      2022-01-06 [1] CRAN (R 4.3.0)
##    class          7.3-22     2023-05-03 [1] CRAN (R 4.3.1)
##    cli            3.6.1      2023-03-23 [1] CRAN (R 4.3.0)
##    codetools      0.2-19     2023-02-01 [1] CRAN (R 4.3.1)
##    colorspace     2.1-0      2023-01-23 [1] CRAN (R 4.3.0)
##    crayon         1.5.2      2022-09-29 [1] CRAN (R 4.3.0)
##    data.table     1.14.8     2023-02-17 [1] CRAN (R 4.3.0)
##    desc           1.4.2      2022-09-08 [1] CRAN (R 4.3.0)
##    devtools       2.4.5      2022-10-11 [1] CRAN (R 4.3.0)
##    dials          1.2.0      2023-04-03 [1] CRAN (R 4.3.0)
##    DiceDesign     1.9        2021-02-13 [1] CRAN (R 4.3.0)
##    digest         0.6.33     2023-07-07 [1] CRAN (R 4.3.0)
##    dplyr        * 1.1.2      2023-04-20 [1] CRAN (R 4.3.0)
##    ellipsis       0.3.2      2021-04-29 [1] CRAN (R 4.3.0)
##    evaluate       0.21       2023-05-05 [1] CRAN (R 4.3.0)
##    fansi          1.0.4      2023-01-22 [1] CRAN (R 4.3.0)
##    farver         2.1.1      2022-07-06 [1] CRAN (R 4.3.0)
##    fastmap        1.1.1      2023-02-24 [1] CRAN (R 4.3.0)
##    forcats      * 1.0.0      2023-01-29 [1] CRAN (R 4.3.0)
##    foreach        1.5.2      2022-02-02 [1] CRAN (R 4.3.0)
##    fs             1.6.3      2023-07-20 [1] CRAN (R 4.3.0)
##    furrr          0.3.1      2022-08-15 [1] CRAN (R 4.3.0)
##    future         1.33.0     2023-07-01 [1] CRAN (R 4.3.0)
##    future.apply   1.11.0     2023-05-21 [1] CRAN (R 4.3.0)
##    generics       0.1.3      2022-07-05 [1] CRAN (R 4.3.0)
##    ggforce        0.4.1      2022-10-04 [1] CRAN (R 4.3.0)
##    ggplot2      * 3.4.3      2023-08-14 [1] CRAN (R 4.3.0)
##    ggpubr         0.6.0      2023-02-10 [1] CRAN (R 4.3.0)
##    ggrepel        0.9.3      2023-02-03 [1] CRAN (R 4.3.0)
##    ggsignif       0.6.4      2022-10-13 [1] CRAN (R 4.3.0)
##    globals        0.16.2     2022-11-21 [1] CRAN (R 4.3.0)
##    glue           1.6.2      2022-02-24 [1] CRAN (R 4.3.0)
##    gower          1.0.1      2022-12-22 [1] CRAN (R 4.3.0)
##    GPfit          1.0-8      2019-02-08 [1] CRAN (R 4.3.0)
##    gprofiler2     0.2.3      2024-02-23 [1] CRAN (R 4.3.1)
##    gtable         0.3.3      2023-03-21 [1] CRAN (R 4.3.0)
##    hardhat        1.3.0      2023-03-30 [1] CRAN (R 4.3.0)
##    highr          0.10       2022-12-22 [1] CRAN (R 4.3.0)
##    hms            1.1.3      2023-03-21 [1] CRAN (R 4.3.0)
##    htmltools      0.5.6      2023-08-10 [1] CRAN (R 4.3.0)
##    htmlwidgets    1.6.2      2023-03-17 [1] CRAN (R 4.3.0)
##    httpuv         1.6.11     2023-05-11 [1] CRAN (R 4.3.0)
##    httr           1.4.7      2023-08-15 [1] CRAN (R 4.3.0)
##    ipred          0.9-14     2023-03-09 [1] CRAN (R 4.3.0)
##    iterators      1.0.14     2022-02-05 [1] CRAN (R 4.3.0)
##    jquerylib      0.1.4      2021-04-26 [1] CRAN (R 4.3.0)
##    jsonlite       1.8.7      2023-06-29 [1] CRAN (R 4.3.0)
##    knitr          1.43       2023-05-25 [1] CRAN (R 4.3.0)
##    labeling       0.4.2      2020-10-20 [1] CRAN (R 4.3.0)
##    later          1.3.1      2023-05-02 [1] CRAN (R 4.3.0)
##    lattice        0.21-8     2023-04-05 [1] CRAN (R 4.3.1)
##    lava           1.7.2.1    2023-02-27 [1] CRAN (R 4.3.0)
##    lazyeval       0.2.2      2019-03-15 [1] CRAN (R 4.3.0)
##    lhs            1.1.6      2022-12-17 [1] CRAN (R 4.3.0)
##    lifecycle      1.0.3      2022-10-07 [1] CRAN (R 4.3.0)
##    listenv        0.9.0      2022-12-16 [1] CRAN (R 4.3.0)
##    lubridate    * 1.9.2      2023-02-10 [1] CRAN (R 4.3.0)
##    magrittr       2.0.3      2022-03-30 [1] CRAN (R 4.3.0)
##    MASS           7.3-60     2023-05-04 [1] CRAN (R 4.3.1)
##    Matrix         1.6-1      2023-08-14 [1] CRAN (R 4.3.0)
##    memoise        2.0.1      2021-11-26 [1] CRAN (R 4.3.0)
##    mime           0.12       2021-09-28 [1] CRAN (R 4.3.0)
##    miniUI         0.1.1.1    2018-05-18 [1] CRAN (R 4.3.0)
##    munsell        0.5.0      2018-06-12 [1] CRAN (R 4.3.0)
##  R myUtilities  * 0.0.0.9000 <NA>       [?] <NA>
##    nnet           7.3-19     2023-05-03 [1] CRAN (R 4.3.1)
##    parallelly     1.36.0     2023-05-26 [1] CRAN (R 4.3.0)
##    parsnip        1.1.1      2023-08-17 [1] CRAN (R 4.3.0)
##    pillar         1.9.0      2023-03-22 [1] CRAN (R 4.3.0)
##    pkgbuild       1.4.2      2023-06-26 [1] CRAN (R 4.3.0)
##    pkgconfig      2.0.3      2019-09-22 [1] CRAN (R 4.3.0)
##    pkgload        1.3.2.1    2023-07-08 [1] CRAN (R 4.3.0)
##    plotly         4.10.4     2024-01-13 [1] CRAN (R 4.3.1)
##    polyclip       1.10-4     2022-10-20 [1] CRAN (R 4.3.0)
##    prettyunits    1.1.1      2020-01-24 [1] CRAN (R 4.3.0)
##    processx       3.8.2      2023-06-30 [1] CRAN (R 4.3.0)
##    prodlim        2023.08.28 2023-08-28 [1] CRAN (R 4.3.0)
##    profvis        0.3.8      2023-05-02 [1] CRAN (R 4.3.0)
##    promises       1.2.1      2023-08-10 [1] CRAN (R 4.3.0)
##    ps             1.7.5      2023-04-18 [1] CRAN (R 4.3.0)
##    purrr        * 1.0.2      2023-08-10 [1] CRAN (R 4.3.0)
##    R6             2.5.1      2021-08-19 [1] CRAN (R 4.3.0)
##    ragg           1.2.5      2023-01-12 [1] CRAN (R 4.3.0)
##    Rcpp           1.0.11     2023-07-06 [1] CRAN (R 4.3.0)
##    readr        * 2.1.4      2023-02-10 [1] CRAN (R 4.3.0)
##    recipes        1.0.8      2023-08-25 [1] CRAN (R 4.3.0)
##    remotes        2.4.2.1    2023-07-18 [1] CRAN (R 4.3.0)
##    rlang          1.1.1      2023-04-28 [1] CRAN (R 4.3.0)
##    rmarkdown      2.24       2023-08-14 [1] CRAN (R 4.3.0)
##    rpart          4.1.19     2022-10-21 [1] CRAN (R 4.3.1)
##    rprojroot      2.0.3      2022-04-02 [1] CRAN (R 4.3.0)
##    rsample        1.2.0      2023-08-23 [1] CRAN (R 4.3.0)
##    rstatix        0.7.2      2023-02-01 [1] CRAN (R 4.3.0)
##    rstudioapi     0.15.0     2023-07-07 [1] CRAN (R 4.3.0)
##    sass           0.4.7      2023-07-15 [1] CRAN (R 4.3.0)
##    scales         1.2.1      2022-08-20 [1] CRAN (R 4.3.0)
##    sessioninfo    1.2.2      2021-12-06 [1] CRAN (R 4.3.0)
##    shiny          1.7.5      2023-08-12 [1] CRAN (R 4.3.0)
##    stringi        1.7.12     2023-01-11 [1] CRAN (R 4.3.0)
##    stringr      * 1.5.0      2022-12-02 [1] CRAN (R 4.3.0)
##    survival       3.5-5      2023-03-12 [1] CRAN (R 4.3.1)
##    svglite        2.1.1      2023-01-10 [1] CRAN (R 4.3.0)
##    systemfonts    1.0.4      2022-02-11 [1] CRAN (R 4.3.0)
##    textshaping    0.3.6      2021-10-13 [1] CRAN (R 4.3.0)
##    tibble       * 3.2.1      2023-03-20 [1] CRAN (R 4.3.0)
##    tidyr        * 1.3.0      2023-01-24 [1] CRAN (R 4.3.0)
##    tidyselect     1.2.0      2022-10-10 [1] CRAN (R 4.3.0)
##    tidyverse    * 2.0.0      2023-02-22 [1] CRAN (R 4.3.0)
##    timechange     0.2.0      2023-01-11 [1] CRAN (R 4.3.0)
##    timeDate       4022.108   2023-01-07 [1] CRAN (R 4.3.0)
##    tune           1.1.2      2023-08-23 [1] CRAN (R 4.3.0)
##    tweenr         2.0.2      2022-09-06 [1] CRAN (R 4.3.0)
##    tzdb           0.4.0      2023-05-12 [1] CRAN (R 4.3.0)
##    urlchecker     1.0.1      2021-11-30 [1] CRAN (R 4.3.0)
##    usethis        2.2.2      2023-07-06 [1] CRAN (R 4.3.0)
##    utf8           1.2.3      2023-01-31 [1] CRAN (R 4.3.0)
##    vctrs          0.6.3      2023-06-14 [1] CRAN (R 4.3.0)
##    viridisLite    0.4.2      2023-05-02 [1] CRAN (R 4.3.0)
##    vroom          1.6.3      2023-04-28 [1] CRAN (R 4.3.0)
##    withr          2.5.0      2022-03-03 [1] CRAN (R 4.3.0)
##    workflows      1.1.3      2023-02-22 [1] CRAN (R 4.3.0)
##    xfun           0.40       2023-08-09 [1] CRAN (R 4.3.0)
##    xtable         1.8-4      2019-04-21 [1] CRAN (R 4.3.0)
##    yaml           2.3.7      2023-01-23 [1] CRAN (R 4.3.0)
##    yardstick      1.2.0      2023-04-21 [1] CRAN (R 4.3.0)
## 
##  [1] /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library
## 
##  R ── Package was removed from disk.
## 
## ──────────────────────────────────────────────────────────────────────────────